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Abstract. Space resolved concentrations of helium He ("^Si) metastable atoms in an 
atmospheric pressure radio-frequency micro-plasma jet were measured using tunable 
diode laser absorption spectroscopy. The spatial profile of metastable atoms in the 
volume between the electrodes was deduced for various electrode gap distances. Density 
profiles reveal the sheath structure and reflect the plasma excitation distribution, 
as well as the dominance of the a-mode discharge. Gap width variations show the 
transition from a normal glow plasma to a pure sheath discharge. In order to analyze 
and verify the experimentally observed profiles of the metastable atoms a 2-dimensional 
simulation model was set up. Applying an appropriate He/N2/02 chemistry model the 
correlation between the metastable profiles and the underlying excitation mechanisms 
was obtained. 
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1. Introduction 

Understanding the energy transfer processes in micro-plasmas is one of the key issues to 
develop new applications and reliable process control. In this context metastable species 
play a decisive role. Due to their long lifetime metastables collide more frequently 
with other particles. The metastable density in micro-discharges is several orders of 
magnitude lower than the density of the ground-state atoms. However, compared to 
most other species the density is significant and the electron collision excitation cross 
sections of some helium levels out of the metastable states exhibit values which are 
several orders of magnitude larger and have much lower thresholds than those for the 
ground state [H [2]. Among these metastable atoms the He (^Si) level plays a decisive 
role, since helium is used as a feed gas in many micro-discharges. This work reports the 
measurement and simulation of He (^Si) atoms in an atmospheric pressure micro-plasma 
jet, providing a simple design by featuring an a-mode RF discharge between two bare 
metallic electrodes. 

Reliable techniques are needed for the systematic investigation of plasma 
characteristics and dynamics. This is essential for the optimization of plasma sources 
and process control. Application of conventional diagnostics, especially invasive ones, 
is often impossible regarding the small dimensions, high operating pressures, and high 
power densities. In this context absorption spectroscopy is a widely used technique 
to measure e.g. absolute concentrations of particles in plasma discharges, since it is 
non- invasive, highly sensitive and provides a sufficient spatial resolution [21 Hj. We have 
applied tunable diode laser absorption spectroscopy (TDLAS) to record the spectral 
profiles of the lowest helium metastable state, deducing absolute densities for various 
discharge conditions |5]. The here used jet configuration is a well analyzed and value 
proven discharge. It was studied already under several experimental and theoretical 
approaches, describing the discharge characteristics in detail P, [71 [8] . The new approach 
of the measurements at hand focuses on the variation of the gap width and determination 
of the sheath/bulk structure and infiuences on the distribution of excited species in the 
discharge. 

To verify the experimental results and to adequately describe the actual excitation 
mechanisms a reliable model is needed. By means of the 2D fiuid dynamic code 
nonPDPSIM ^ we perform simulations which take into account both the discharge 
dynamics between the electrodes as well as the infiuence of the lateral gas fiow. 

2. The atmospheric pressure micro-plasma jet 

The atmospheric pressure micro-plasma jet is a capacitively coupled, non-thermal glow- 
discharge plasma at high pressures. The design concept of this discharge is based on 
the plasma jet introduced by Selwyn et al. in 1998 [ID] and advanced by Schulz-von der 
Gathen et al. [6] termed yuAPP J. The feed gas fiows between two closely spaced stainless 
steel electrodes driven at 13.56 MHz radio-frequency in a parallel plate configuration 
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(Figure [T]). Electrodes, plasma volume and effluent are enclosed by quartz windows, 
giving direct optical access to the plasma itself and the effluent volume behind the 
electrodes. The discharge uses helium as feed gas with typical gas velocities around 
100 ms~^. The electric field between the electrodes causes a breakdown in the gas 
and produces a plasma with electron temperature and density of about 1 to 2eV and 
10^°cm~^, respectively [71 |8]. Atoms and molecules in the feed gas become excited, 
dissociated or ionized by electron impacts. Since the electrons are not in thermal 
equilibrium with the ions and neutrals, the gas temperature remains a few tens Kelvin 
above room temperature 

The presented jet configuration features a dielectric extension of the gas channel to 
assure controlled gas flow in the effluent behind the plasma. Electrodes and dielectric 
extensions are 40 mm and 50 mm in length, respectively. The distance between the two 
windows is fixed to 1mm while the electrode gap width is variable between 0.2 mm and 
30 mm. 

2.1. TDLAS set-up 

The small dimensions of micro-discharges and their operation at atmospheric pressures 
are a challenge for optical diagnostics, since high sensitivity and high spatial resolution 
are required. For the absorption spectroscopic measurements a standard TDLAS set-up 
was used. The absorption profile was recorded by scanning the laser frequency across the 
absorption line. Figure [2] shows a sketch of the experimental setup. The laser beam from 
the DL passed through two beam splitters. A part of the beam was guided to a Fabry- 
Perot interferometer (1 GHz free spectral range), a second part through a low pressure 
reference cell to perform the calibration of the laser frequency. The part of the beam 
transmitted through the first beam splitter was attenuated by neutral density filters with 
an optical density in the order of 3, and focused into the discharge with a beam power 
of less than 2 /iW at 100 /im spot size, to avoid any saturation effects. After passing 
the discharge the beam was guided through a set of apertures and filters to suppress 
the emission from the plasma by reducing the collection angle and blocking wavelengths 
different than the observed transitions. The transmitted beam intensity was measured 
by photodiodes with on-chip transimpedance amplifier. For an effective measurement of 
the absorption signal across the jet axes, the discharge casing was mounted on a small 
movable stage featuring three electronically controlled stepping motors to adjust, with 
high precision, the discharge cell position in all spatial dimensions. This set-up allows 
the positioning of the jet with an accuracy of about 5 /xm and automated a;2;-mapping 
of the complete plasma volume (In y direction all measurements presented in this paper 
are space averaged). 

Since the absorption rate of the laser light by metastable atoms is very low, in the 
order of 10~^ after 1 mm absorption length, lock-in technique was used to measure the 
changes in signal intensity. Applying lock-in technique requires the pulsing of the signal 
to be measured. This was realized by pulsing the RF-power coupled into the system. 
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which consequently leads to a pulsing of the metastable density in the discharge. The 
pulse frequency was chosen to 4 kHz and a duty cycle of 50 %. 

Absolute metastable densities were derived from the transmittance I/Iq and the 
Beer-Lambert law. Therefore, the four signals 

L(z/) — Plasma and Laser on, 

Lq^u) — Plasma off. Laser on, 

P(z/) — Plasma on. Laser off, 

-B(z/) — Plasma and Laser off (background), 

have to be acquired to calculate the transmittance spectra and correlate them with the 
plasma properties by 

Loiu) - B{v) ■ ^ ' 

and /o(^) are the intensities of transmitted radiation with and without the presence 
of absorbing species, k{i') is the absorption coefficient and / the path length through 
the absorbing medium [12]. The absorption coefficient is connected to the population 
density of metastable atoms by 

m = -^^hkN,F{v), (2) 
47reo cmg 

where fit is the oscillator strength of the line, Ni the density of the lower level, and 
F{v) a normalized function F{v)dv = 1) representing the absorption lineshape. All 
other terms have their usual definitions. The absolute metastable density can then be 
given by 

where 5* is the area under the absorption curve that provides the line-averaged density 
of the absorbing species. 



oo 



2.2. Simulation set-up 

The 2D simulations are performed on a non-structured numerical grid with the fluid 
dynamic code nonPDPSIM which has been originally designed and realized by Kushner 
and co-workers. A detailed description of the code and a number of successful 
applications can be found in [HI HSj [TH [I5]. Basic simulations of a radio- frequency 
driven He/02 plasma jet and its effluent accounting for the plasma dynamics and the 
lateral gas flow can be found in [16]. Here, we just briefly discuss the implemented 
equations and the underlying physics. 

The code nonPDPSIM simulates the dynamics of weakly ionized plasmas in the 
regime of medium to high pressure. It takes into account the physics and chemistry 
of charged particles - electrons with mass and charge — e, ions with mass nij and 
charge qj - and of the excited as well as the ground state neutrals (mass mj). For all 
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species j, the continuity equations (particle balances) are simultaneously solved, where 
Tj is the particle flux density and Sj is the source and loss term, respectively: 

^ = -V-f, + 5,. (4) 

The fluxes are calculated from the momentum balances in the drift- diffusion 
approximation evaluated in the local center-of-mass system. Dj and nj are the diffusion 
constant and the mobility (if applicable) of species j. Further, E is the electrical field, 
and V is the mass-averaged advective velocity of the medium: 

= n^v - DjVn^ + ^^fi^rijE. (5) 

For the electron fluid, additionally an energy balance equation is solved taking into 
account Ohmic heating and the energy losses due to elastic and inelastic interaction 
with the neutrals and ions as well as heat conduction, 

^ (^^neTe^ = i-E-S7 ■ (^-KeVTe + <^ " "-e X! ^(^ihrii. (6) 

To capture the non-Maxwellian behavior of the electrons, all electronic transport 
coefficients (the mobihty /Xe, the diffusion constant De, and the thermal conductivity 
Kg) as well as the electronic rate coefficients in eqs. ^ and (jSD are calculated by the 
local mean energy method: A zero- dimensional Boltzmann equation for the electron 
energy distribution /(e) and the transport and rate coefficients is solved for the locally 
applicable gas composition and various values of the electrical field. The tabulated 
data are then consulted in dependence of the fluid dynamically calculated electron 
temperature T^. 

The plasma equations are coupled to a modifled version of the compressible Navier- 
Stokes equations which are solved for the gas density p, the mean velocity and the gas 
temperature T . The contributions to the energy equation from Joule heating include 
only ion contributions; the heat transfer from the electrons is included as a collisional 
change in the enthalpy. The scalar pressure p is given by the ideal gas law. 

Finally, the potential $ is calculated from Poisson's equation. (The code works in 
the electrostatic approximation such that E = — V$.) The charge density on its right 
hand side stems from the charged particles in the plasma domain and from the bound 
charges ps at the surfaces. The coefficient e = £o^r represents the permittivity of the 
medium: 

-V-eV^ = Y.qjnj + p,. (7) 

j 

The surface charges are governed by a separate balance equation, where a is the 
conductivity of the solid materials and the subscript s indicates evaluation on the 
surface: 



dps 
dt 



Eg,(-V-f, + 5,)-V-(a(-V$)) 



(8) 
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The dynamical equations are complemented by an appropriate set of boundary 
conditions. Electrically, the walls are either powered or grounded. With respect to the 
particle flow, they are either solid, or represent inlets or outlets: The flow is specified to a 
given flux, while the outlet flow is adjusted to maintain the pressure. Finally, it is worth 
mentioning that the actual implementation of the equations poses some difficulties due 
to the vast differences in the time scales of the dynamics of the plasma and the neutrals. 
These difficulties are overcome by the methods of time-slicing and subcycling. 

The described code is employed to simulate the /zAPP J depicted in Figure [1] with 
varying electrode gap sizes of 1.8 mm, 1.0 mm and 0.3 mm. We choose these values to 
characterize the three distinct regimes of plasma bulk-sheath ratios. The simulation 
resolves the two Cartesian dimensions x and z. In direction y, translational invariance 
is assumed. 

By leaks in the gas supply system as well as intrusion of air from the exit nozzle 
of the jet, the amount of nitrogen and oxygen in the feed gas is significant and has to 
be taken into account for the modeling approach. We mimic the impurities of the 
experimental setting according to [17] by adding 0.016 % of nitrogen and 0.004 % 
oxygen in the chemical model (Oxygen species contribute to the electronegativity of 
the plasma while nitrogen has an impact on the quenching of the helium metastable 
density.). We consider the following species: Ground state neutrals O2, O, O3, and 
He, 02(1'), representing the first four vibrational levels of O2, the electronically excited 
states 02(^Ag), 02(^S+), 0{^D), 0(^5), and He* = He(3Si), He^ = He2(^E+), positive 
ions O2, O"*', He"*" and He^, negative ions O^, 0~, and O3, and electrons. The 
flow rate of the gas mixture injected in the jet is adjusted to achieve advective gas 
velocities of about 100 ms~^. The outlet is controlled to maintain a constant pressure. 
Finally, the secondary electron emission coefficient of the electrode surfaces depends on 
the ion species and varies from 0.26 for He^ to 0.06 for 

3. Results and discussion 

3.1. Spatial distribution of He metastable atoms 

Figure [3] shows 2D-maps of the He (^Si) metastable density in the discharge volume. The 
upper map shows experimental results whereas the bottom map shows simulation data. 
The horizontal and vertical axis span the exact area between the electrodes. The map 
of the experimental data, which covers 2,000 reading points (40 vertical x 50 horizontal) 
of the absorption signal, in the plasma volume is in good qualitative and quantitative 
agreement with the simulation results. It should be mentioned that the steps occurring 
in the measured map at certain positions are due to mechanical limitations in the setup, 
that made it necessary to record that map in four successive intervals. Furthermore 
the measurement time for the map was several days, explaining why the map shows 
heavy fluctuations, whereas the actual metastable distribution in the discharge is more 
smooth, as seen in the measurements presented later in this article. 
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A variety of effects determines tlie metastable distribution in the horizontal and the 
vertical axis. In horizontal direction the density profile is governed mainly by two effects. 
One is the alignment of the electrodes, since the electric field distribution strongly 
influences the electron energy. In the present experimental case the electrode gap is 
marginally larger on the rear side of the jet, resulting in a slight decrease of densities from 
left to right. A second impact is given by impurities entering the jet through the exit 
nozzle, therewith quenching the metastable atoms. Since the jet is running permanently 
in contact with the ambient atmosphere the intrusion of nitrogen and oxygen into the 
plasma channel is significant, and decreases the metastable density especially in the first 
millimeters from the nozzle. The small fluctuations in the horizontal profile (which are 
not resolved by the simulation model) can be attributed to turbulences, imperfections 
in the surface of the electrodes or singularities in the walls of the gas channel. 

The vertical metastable profile is less determined by these environmental factors 
but by the electron density and temperature distribution. The observed structure is 
consistent with the sheath structure in an RF discharge. Directly in front of the 
electrodes, the metastable density is low since the electron density is too low for an 
efficient excitation and ionization of the ground state atoms (see also Figure HJ below). 
Metastables reach their highest densities some 10 fim away from the electrode at the 
plasma-sheath interface, where in the negative glow area most excitation and ionization 
processes occur. Here in the sheath/bulk interface the electron energy is highest. In 
the bulk, the electric field strength and therewith the electron temperature is low (see 
also Figure [5l below), yielding to a lower metastable production rate. The asymmetry 
in the vertical profile of the simulated data can be attributed to stochastic effects in the 
simulation since it is not averaged over many RF-cycles due to runtime considerations, 
while the measurements are averaged over several 10^ cycles. 

3.2. Gap width variations 

Due to the variability of the electrode gap the micro-plasma jet offers the possibility 
to continuously tune the sheath/bulk ratio and therewith the vertical excitation profile. 
Figure H] shows the He metastable profile for gap widths between 2.7 and 0.25 mm. 
These values represent the limits where the pure helium discharge can be operated in 
a-mode. With decreasing gap width the breakdown voltage decreases while approaching 
the Paschen minimum. Although the outer power coupling to the discharge for each gap 
width is different, one can clearly observe a systematic increase of the metastable density 
in the maximum with smaller gap width. A quantitative description of this increase is 
difficult, since the gap width correlates with a variety of discharge parameters, like the 
impedance, the power coupling or the gas velocity, that all influence the species' density. 

The graphs show a steep rise of the density in front of the electrode surface with a 
maximum in the sheath/bulk interface, while in the bulk the density is low. The den- 
sity distribution suggests that diffusion processes play a minor role for the metastable 
profile. At atmospheric pressure the mean free path of the electrons is short. High 
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energetic electrons are lost in the sheath and can not contribute to the excitation in the 
plasma bulk. Although the metastable lifetime is comparatively long, in the order of 
1 fis, the metastable diffusion to the bulk can be neglected likewise [T7]. De-excitation 
of these species is dominated by three-body collision with ground state atoms as well 
as quenching reaction with impurities, like molecular nitrogen and oxygen. A detailed 
analysis of metastable lifetimes in this discharge is given in [17]. For gap widths larger 
than 1mm the two excitation areas are clearly separated by the plasma bulk. Lower- 
ing the electrodes distance leads to the coalescence of the two sheath regions, and an 
excitation profile where highest densities are in the center of the gap. For electrode 
distance of just a few hundred fim the plasma transforms to a pure "sheath discharge" . 
These experimental findings are supported by the results from the numerical simula- 
tion. Figure |5] (top) shows the metastable profile for selected gap widths, the related 
electron densities and electron temperatures are given in Figure [5] (bottom). Due to 
the admixture of oxygen into the feed gas the electronegativity of the plasma has to be 
taken into account. The electron temperature profile for a gap widths of 1.8 mm rep- 
resents the excitation pattern of the metastable atoms. The temperatures are highest 
in the sheath/bulk interface, with a maximum of about 3eV. Lowering the gap widths 
results, like in the experiment, in a convergence of the two sheaths. For a low gap 
widths of 0.3 mm the electron density increases significantly, while temperatures reach 
values of almost 4eV. This behavior suggests that the high energy tail of the EEDF 
is getting more pronounced. Since the high electron energies enhance the production 
of metastable atoms, the simulation again supports the qualitative observation of very 
high metastable densities for low gap widths. 

For metastable density measurements no error bars are shown, since statistical as 
well as systematic errors are negligible. For the presented data they are in the order 
of lO^cm"^, while the density values are in the order of 10^^ cm~^ to lO^^cm"^. Due 
to the use of lock-in technique a sensitive measurement of low densities was possible 
with a high signal to noise ratio. A minor systematic error is caused by the unknown 
metastable profile in y-axis, and the implications on the absorption length. But since the 
metastable production and loss processes take place on a very local scale at atmospheric 
pressure, the profile in y-axis is expected to be homogeneous with the exception of a 
tiny volume in front of the glass windows that confines the plasma channel and acts as 
a sink for metastable atoms. 

4. Summary 

We have shown that the experimental results for metastable distributions in micro- 
plasma jet discharges can be reproduced qualitatively an quantitatively by a 2D 
numerical model, therewith revealing the underlying excitation mechanisms of the 
discharge. By varying the gap width it was proven that the simulation model reliably 
mimic the experimental results for a large span of sheath bulk ratios. This offers now 
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new possibilities in the analysis of extreme discharge conditions, like 7-mode structures 
and very small gap widths showing effects like no quasi-neutrality all over the discharge 
volume. 
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Figure 1. Sketch of the micro-plasma jet discharge, which is also the simulation 
domain. Shown is a 2-dimensional cross-section through plane that is spanned by the 
two electrodes. 
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Figure 2. TDLAS setup. 
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Figure 3. Top: 2-dimensional map of the measured metastable density in the 
discharge volume for the He — > 2^Pj 2 transition. Powered electrode at the 
top, grounded electrode at the bottom. Electrode gap width was 1.8 mm. Densities 
are given in cm~'^. Bottom: Simulation results for the spatial distribution of the He 
metastables. 
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Figure 4. Measured vertical metastable profile for gap widths between 2.7 mm and 
0.25 mm. 
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Figure 5. Top: Simulation results of characteristic metastable profiles for various 
gap widths. Bottom: Simulation results of electron density and electron temperature 
profiles for various gap widths. 



